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Abstract. We investigate the X-ray properties of the powerful 
radio galaxy Hercules A (3C 348) using ROSAT HRI, PSPC 
and ASCA observations. The ASCA data are well fit by a ther- 
mal plasma model with a temperature of about 4.3 keV and 
abundances 0.4 solar The HRI observation clearly reveals ex- 
tended and elongated X-ray emission. For radii greater than 10 
arcsec, the surface brightness profile perpendicular to the elon- 
gation is well fit by an isothermal /3-model with /? = 0.63 and a 
core radius of « 120 kpc. We derive a total mass of the putative 
cluster of 8.4 x 10^^ Mq and a gas mass fraction of about 18 
per cent. The inner part of the surface brightness profile cleaky 
reveals the presence of a point source, which contributes «8 per 
cent to the total flux. The 0.1-2.4 keV luminosity of the point- 
like and the extended emission is 3.4 x 10^^ and 4.3 x 10'*''^ 
erg s~^, respectively. After subtracting the cluster X-ray emis- 
sion from the HRI image, residual structures are visible, which 
partly coincide with the radio jet and lobes. This indicates an 
interaction of the radio jet with the intracluster medium. 
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1. Introduction 

The giant elliptical radio galaxy Hercules A (3C 348, 1648H-05, 
z = 0.154) is one of the most prominent AGN in the sky with 
a total radio flux density of 351 Jy at 178 MHz, which makes 
it the fourth brightest extragalactic radio source. It was one of 
the first optically identified extragalactic radio sources (Bolton 
1948) and has been studied extensively in the radio and optical 
since then. 

A high resolution radio map of Hercules A at 5GHz (Dreher 
& Feigelson 1984) shows giant radio lobes and prominent jets 
with spectacular (and still unexplained) ring-like structures. 
The classification of Hercules A in terms of Fanaroff & Riley 
(FR) type is ambiguous. Despite its extremely high radio power 
(log Ptot,5GHz = 27.19 W/Hz), which is typical for FR II radio 
galaxies, the radio morphology argues for a FR I classification 



Correspondence to: J. Siebert, jos@mpe.mpg.de 



(no hot spots visible in the radio lobes, prominent two-sided 
jets). 

Broadband optical imaging revealed two concentrations of 
Ught. Thus, Hercules A was originally thought to be a cD 
galaxy with a foreground star superimposed on it (Greenstein 
1962), although already Minkowski (1957) classified it as a 
double galaxy. Recent photometric work by Sadun & Hayes 
(1993) supports the idea of a double nucleus with a separation 
of ~ 4", thus indicating merging activity. 

The environment of Hercules A has been investigated by 
several authors (Greenstein 1962, Yates et al. 1989, Allington- 
Smith et al. 1993), but whether or not Hercules A is associated 
with a cluster is not yet unambiguously decided. Whereas Yates 
et al. (1989) claim no particularly rich environment around 
Hercules A, Allington-Smith et al. (1993) find an excess num- 
ber of galaxies within 500 kpc of Hercules A. It has to be noted, 
however, that all studies were based on statistical methods, i.e. 
number counts, and there is no spectroscopic confirmation of a 
cluster of galaxies up to now. 

In X-rays Hercules A was observed for the first time with 
the Einstein observatory. Dreher & Feigelson (1984) note ex- 
tended emission on scales of > Mpc and a luminosity of 3.4 x 
10^^ erg s~^ in the 0.2^ keV energy band. However, a detailed 
analysis of the X-ray data has never been published. 

In this paper we present a thorough analysis of the X-ray 
properties of Hercules A using archival ROSAT data and new 
ASCA observations. We will show that the X-ray emission is 
extended and predominently thermal and that there is interac- 
tion between the radio jet and the X-ray emitting gas. The pa- 
per is structured as follows: In Sect.2 we analyse the ASCA 
and ROSAT PSPC data and determine the spectral properties 
of Hercules A. In Sect. 3 we investigate the spatial structure 
of the X-ray emission using a ROSAT HRI observation and 
we compare it to the VLA radio map. A discussion of the re- 
sults is given in Sect.4 and our conclusions are presented in 
Sect.5. Spatial scales and luminosities are calculated assuming 
i/o = 50 km s"^ Mpc"\ go = 0.5 and A = throughout this 
paper At the distance of Hercules A one arcmin corresponds to 
w 209 kpc. 

2. Spectral analysis 
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2.1. TheASCA and ROSAT PSPC observations 



The ASCA observation was performed on August 12 and Au- 
gust 13 1998 in 1-CCD faint mode. The data were analysed 
using FTOOLS 4.1. The recommended standard screening crite- 
ria were applied. 

In addition, data taken within 60 seconds after passage 
of the day-night-terminator and the South Atlantic Anomaly 
(SAA) were not considered in the analysis. Periods of high 
background were manually excluded from the data by checking 
the light curve of the observation. The resulting effective expo- 
sures were 35. 1 and 34.4 ksec for G1S2 and G1S3 and 33.4 and 
32.9 ksec for SISO and SlSl, respectively. 

Source counts were extracted from a circular region cen- 
tered on the target with a radius of 6' for the GIS and 4' for 
the SIS. We used the local background determined from the 
observation in the analysis for both detectors. In particular, the 
GIS background was estimated from a source free region at the 
same off-axis angle as the source and with the same size as the 
source extraction region. In total, ^4000 source photons were 
detected in each of the SIS detectors. 

All spectra were rebinned to have at least 20 photons 
in each energy channel. This allows the use of the 
technique to obtain the best fit values for the model spec- 
tra. We used the latest GIS redistribution matrices avail- 
able (V4_0) from the calibration database and created the 
SIS response matrices for our observation using SISRMG, 
which applies the latest charge transfer inefficiency (CTI) table 
(sisph2pi_110397.fits). The ancillary response files for 
aU four detectors were generated using the ASCAARF program. 

Spectra were fitted in the energy range 0.8 to 9 keV for 
both GIS. For SISO and SISl we used 0.8 to 8 keV. The upper 
energy boundaries are given by the maximum energy at which 
the source was detected in each instrument. The lower energy 
boundaries result from the calibration uncertainties of the de- 
tectors. In particular, SIS 1 recently shows systematic residuals 
below 0.6 keV (Dotani et al. |1997[). Also, since the CCD tem- 
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perature of both SIS detectors was rather high during the ob- 
servation, we decided to ignore all data below 0.8 keV. An in- 
creasing RDD (Residual Dark Distribution), which is caused by 
radiation damage, has been reported for both CCDs. This effect 
cannot be corrected for with currently available software. How- 
ever, RDD degradation should be negligible for 1-CCD obser- 
vations and affects mostly the lowest energy channels, which 



are excluded from our analysis anyway (Dotani et al. 1997). 

The ROSAT PSPC data (ROR 70161 1) were obtained from 
the public data archive at MPE. Hercules A was observed with 
the ROSAT PSPC between August 19 and August 3 1 1993 with 
an effective exposure of 8058 sec. The data were prepared for 
spectral analysis using standard commands within the EXSAS 
environment. In short, photons were extracted from a circu- 
lar region with a radius of 5 arcmin centered on the peak of 
the X-ray emission. The background was determined from a 
source-free annulus with inner radius 6 arcmin and outer ra- 
dius 9 arcmin. Only pulse height channels 12 to 240 were used 
for spectral analysis, because of calibration uncertainties in the 
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Fig. 1. Ratio of the ASCA data from all four detectors to a sim- 
ple power law model (top panel) and the best-fit thermal model 
including a five per cent power law contribution from the AGN 
(bottom panel). The data are rebinned to match the energy res- 
olution of the instruments. 

lowest energy channels. The data were rebinned to get a signal- 
to-noise ratio of at least five in each energy bin. Finally, vi- 
gnetting and dead time correction were applied to the binned 
dataset. 

2.2. Results 

We first investigated the data from each detector individually 
and the results agreed within the errors. In the following we 
therefore only cite the best-fit model parameters of a simulta- 
neous fit to all four instruments. 

The ratio of the ASCA data to a simple power law model 
with Galactic absorption (iVn — 6.40 x 10^° cm^^) is shown 
in the upper panel of Fig. |l|. Clearly, there are systematic resid- 
uals, in particular above 3 keV. For completeness, we give the 
best-fit photon index: T = 1.9ltg:o3 (x^ = 645.4 (586 d.o.f.)). 

Next we tried a thermal emission model, namely the 
Mekal model within XSPEC. The fit improves dramatically 
(Ax^ = —61.7 for one additional parameter). The best-fit pa- 



rameters are as follows: kT = 4.62 



+0.27 
-0.25 



keV, Z = 0.38 ± 



O.1GZ0 and A = (5.0±G.2) x lO^^ photons cm^^ g-i keV^i. 
Since the spatial analysis indicates a small contribution from 
the central AGN (see Sect.3), we finally included a power law 
component and fixed its normalization to 10 per cent of the 
thermal emission. The fit does not further improve (Ax^ = 
+1.1 for one additional parameter), but the model is well con- 
sistent with the data (Xrcd ^ 1-00). The best-fit parameters of 



this model are kT = 4.25 



-1.00 
-0.66 



keV, Z = 0.44 ± O.UZr: 
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r = 1.7i;^;^o and Ah = ^^-O^to.lr x lO'^ photons cm-^ 
keV~^. The corresponding data to model ratio is shown in 
the lower panel of Fig. |l]. 

If the above model is applied to the PSPC data with all pa- 
rameters fixed to their best-fit values, we only get a marginally 
acceptable description of the PSPC spectrum (Xrod = ^-"^^ 
(37 d.o.f.))- However, the fit significantly improves for lower 
temperatures (kT — 2.44lQ'gQ keV; all other parameters fixed 
to the ASCA values). This might indicate that a single tem- 
perature thermal plasma is not a good description for the in- 
tracluster medium around Hercules A. In the case of a more 
complex temperature structure, the relative contribution of the 
lower temperature plasma in the ROSAT band is higher than in 
the ASCA band. Therefore the PSPC observation is more sen- 
sitive to the lower temperature component, whereas the ASCA 
measurement is dominated by the higher temperature compo- 
nent. However, the fact that the best-fit temperature determined 
from the PSPC data is close to the upper end of the PSPC en- 
ergy range might also indicate that instrumental effects and 
maybe even residual ROSAT-ASCA cross-calibration uncer- 
tainties contribute to the difference in the obtained plasma tem- 
peratures from the two instruments. 

We further investigated if there is a radial temperature gra- 
dient, but found none. Fitting the PSPC data from the inner 45 
arcsec and from an annulus with inner radius 45 arcsec and 
outer radius 3 arcmin seperately, we do not find significantly 
different temperatures (/cTjn — 2.17^^'^^ keV compared to 
fcTout = 2.68lo 72 keV). Also the (lower resolution) ASCA 
data do not indicate any temperature gradient. 

3. Spatial analysis 

3.1. The HRI observation 

Hercules A was observed with the ROSAT High resolution Im- 
ager (HRI) between August 28 and September 1 1 1996 ( ROR 
702755). The effective exposure of this observation was 21.52 
ksec. A second, much shorter (1.2 ksec) observation is not con- 
sidered in this analysis. The total count rate for Hercules A is 
0.077 ± 0.010 cts/s. 

Only HRI pulse-height channels 2 to 8 are used in the spa- 
tial analysis in order to increase signal-to-noise, as channels 9 
to 15 mostly contain instrumental background (cf. David et al. 
1997). A contour plot of the total X-ray emission is shown in 
Fig. H The photons were binned in 2 x 2 arcsec pixels and 
subsequently smoothed with a Gaussian with cr = 6 arcsec. 

The X-ray emission is clearly extended and elongated. As 
will be shown later, the direction of the elongation is close to 
the radio axis of Hercules A. We are confident that the peculiar 
shape of the X-ray emission is indeed intrinsic and not an in- 
strumental effect, because there is a second X-ray source about 
4 arcmin to the south-east of Hercules A (R.A. = 16h5 lm22.3s, 
Dec = H-04d58m23s), which is point-like. It is associated with 
a stellar object of my = 12.77 on the digitized POSS-I plate. 
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Fig. 2. X-ray contours of Hercules A. The contours correspond 
to 5,8,12,16,20,30 and 50 a above background. The lowest 
and highest contours denote intensity levels of 9.4 x 10"'^ and 
61.8 X 10^^ cts s^^ arcmin^^ respectively. 

3.2. Surface brightness profile 

To determine the physical properties of the extended X-ray 
emission we fitted a /3-model (e.g. Cavaliere & Fusco-Femiano 
1976; Jones & Forman 1984) of the form 

/ ^2 X -3/3+1/2 

to the surface brightness profile of the HRI source. However, 
we did not use an azimuthally averaged radial profile, but ex- 
tracted photons from a cone with an opening angle of 90 de- 
grees oriented perpendicular to the elongation of the X-ray 
emission. Since we intend to subtract the pure cluster emission 
to see if the residuals are correlated with the radio structure, the 
above procedure avoids subtracting too much of the extended 
X-ray emission which originally might be associated with the 
radio source. 

The surface brightness profile, shown in Fig. ||, peaks at 
the position of the AGN. In the fitting procedure we there- 
fore only used data points for radial distances greater than 10 
arcsec to avoid contamination by AGN emission. The best-fit 
values for the /3-model are So = 647.1^200 5 cts/arcmin^, 
/3 = 0.631qo3 and Vc — 34.6l;3'3 arcsec (corresponding to 
120.4j^![5 5 kpc). To account for the central AGN emission we 
included the Point-Spread-Function (PSF) model (David et al. 
1997) for the HRI in the fit. The theoretical PSF-model was 
convolved with an additional Gaussian (cr — 1.5 arcsec) to 
take into account the known smearing of the PSF by residual 
wobble motion, which varies between different observations 
(Morse 1994). For the normalisation of the PSF model we get 
3973.8l^QgQQ cts arcmin^^. By integrating the two profiles we 
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Fig. 3. The surface brightness profile for Hercules A. The best- 
fit /3-model is indicated by the long-dashed line, whereas the 
background level and the the PSF model for the central point 
source are given by the dotted and the dot-dashed lines respec- 
tively. 



calculate the contribution of the AGN to the total X-ray emis- 
sion to be about 8 per cent. Assuming a power law spectrum 
with r = 1.7 and Galactic absorption we get a 0. 1-2.4 keV flux 
from the point source of fx ~ 3.3 x lO^^'^ erg s^^ cm^^ which 
corresponds to a rest frame 0.1-2.4 luminosity of 3.4 x lO"'"^ 
erg s~^. Hercules A has been noted to be over-luminous in X- 
rays with respect to its optical continuum and line luminosities 
(Siebert et al. 1996; Brinkmann et al. 1995). However, these 
investigations were based on the ROSAT All-Sky Survey data 
for Hercules A and no separation of extended and AGN emis- 
sion could be done. The much lower AGN luminosity obtained 
from our spatial analysis places Hercules A with the bulk of 
radio galaxies investigated in the above mentioned studies. The 
flux and luminosity for the extended emission component are 
/x « 4.1 X 10^12 erg cm^^ and ~ 4.3 x lO'^'* erg s^^ 
(assuming a thermal Bremsstrahlung spectrum of kT — 4.3 
keV; see below). The extended X-ray luminosity is typical for 
clusters of galaxies with Abell richness class two or higher 
(Ebeling 1993). This would also be consistent with the number 
counts of AUington-Smith et al. (1993) to determine the galaxy 
density around Hercules A. In addition, the observed luminos- 
ity and temperature for Hercules A are consistent with the well 
known luminosity-temperature relation for clusters of galaxies 
(David et al. 1993). 

To investigate the residual X-ray emission from the AGN 
we subtracted the best-fit one dimensional /^-model given 
above from the binned and smoothed HRI image. The resulting 
X-ray contours overlayed on an optical image from the corre- 
sponding POSS-I plate are shown in Fig. ^ 

The residual X-ray emission is dominated by the point 
source which is most likely associated with the active nucleus 
of Hercules A. Several other sources contribute to the remain- 
ing X-ray emission. A significant fraction of the diffuse com- 



Fig. 4. Residual X-ray contours after subtracting the best-fit /3- 
model from the HRI image overlayed on an optical image of 
the field from the digitized POSS-I plates. The contours corre- 
spond to 3,5,8,12,16,24 and 32 a above background. The low- 
est and highest contour denote intensity levels of 7.0 x 10~^ 
and 40.3 x lO^'' cts s^^ arcmin^^ respectively. 



ponent might still be unsubtracted cluster emission, since we 
assumed a circular symmetry in the cluster subtraction pro- 
cedure. Part of the emission might also come from additional 
point sources, for example the bright (my — 11.97) stellar ob- 
ject to the northwest of Hercules A. Other features, as for exam- 
ple the one-sided jet-like emission and the diffuse component 
about 1 .5 arcmin to the west of Hercules A, correlate with the 
radio emission. They will be discussed in detail in Sect. 3.4. 

3.3. Physical parameters of the Cluster emission 

By deprojection of the surface brightness profile one can derive 
the corresponding density profile: 



(1) 



n{r) = no 1 + 



The central density uq is then given by (e.g. Henry et al. 1993): 



no = 1-2 X lO^^cm"^ 



1/2 



r(3/3 - 1/2) 

r(3/3) 



7 0.7, 



E2 

kT 



-7 0.7, 



El 

kT 



-1/2 



Iq is the central intensity 5*0 converted to erg cm^^ s^^ sr~^, 
Tc the core radius in cm and kT is given in keV. F and 7 are 
the complete and incomplete Gamma functions. Ei and E2 de- 
note the observed energy range, i.e. 0.1 and 2.4 keV in our 
case. For the conversion from 5*0 to Iq we assumed a thermal 



Siebert et al.: Hercules A 



5 



bremsstrahlung spectrum with kT = 4.3 keV and Galactic ab- 
sorption. Using the best-fit parameters from our /3-model we 
obtain uq = 9.1 x 10~^ cm^''. By integrating Eq. (1), we can 
now calculate the total gas mass within a given radius. In the 
case of Hercules A we get 1.6 x 10^"^ Mq within six times the 
core radius, i.e. about 600 kpc. 

Assuming hydrostatic equilibrium for the intracluster 
medium of Hercules A and spherical symmetry, we derive the 
total gravitating mass as a function of radius from the tempera- 
ture and gas density profiles by applying the hydrostatic equa- 
tion: 



Mtot(r) 



dp{r) 



kTg{r)r 

mnj^G \Tg{r) dr p{r) dr 



(2) 



As discussed in Bohringer et al. (1998), deviations from hydro- 
static equilibrium as well as moderate ellipticities do not have a 
large effect on the mass determination. With the additional as- 
sumption of an isothermal intracluster medium, Eq. (2) reduces 
to 



Mtot(r) 



3(3kTg 



G^imurl 1 + ^3 



(3) 



For the total gravitating mass within 600 kpc we get 8.4 x 10^^ 
Mq. Hence, the gas mass fraction is w 18 per cent at this ra- 
dius. 



3.4. Radio/X-ray interaction 

Fig. ^ shows the radio contours from a 5GHz VLA observation 
(Dreher & Feigelson 1984), overlayed onto the HRI image af- 
ter subtraction of the cluster emission according to the best-fit 
/3-model. Several features of the residual X-ray emission coin- 
cide spatially with the radio emission: (a) faint knots of X-ray 
emission are surrounding the outer edges of both radio lobes; 
(b) there is a bright patch of X-ray emission close to the head 
of the western radio lobe. It is only slightly displaced from a 
faint radio feature, which might be interpreted as the head of 
the advancing radio jet; (c) the orientation of the jet-like fea- 
ture in the X-ray emission close to the central AGN is aligned 
with the radio jet; (d) the brightest spot in the eastern radio jet 
is bracketed by enhancements in the X-ray emission. We note 
that all features are robust in the sense that they also show up 
when we subtract the cluster emission determined from the az- 
imuthally averaged radial profile. 

The close association of X-ray and radio features strongly 
argues for interaction of the relativistic gas in the radio jet and 
the thermal gas of the intracluster medium. An X-ray cavity co- 
inciding with the radio lobes has been noted previously for the 
radio galaxies NGC 1275 (Bohringer et al. 1993) and Cygnus 
A (Carilli et al. 1994). In the latter case also an enhancement 
of the X-ray emission close to the edges of the radio lobes 
was reported. The general physical scenario for these effects 
(cf. Carilli et al. 1994 and Clarke et al. 1997 for a more 
detailed treatment) is such that the advancing and expand- 
ing jet is ploughing its way through the intracluster medium 
(ICM) thereby expelling the thermal gas from its interior A 



thin sheath of dense, shocked material develops in the vicin- 
ity of the advancing head and the expanding tails of the ra- 
dio lobes. When our line of sight is tangential to the edges of 
the radio lobes, these density enhancements become visible in 
X-rays, because thermal X-ray emission is proportional to the 
square of the gas density. For lines of sight through the center of 
the radio lobes the X-ray enhancements are probably balanced 
by an X-ray cavity within the radio lobes due to the expelled 
gas. The faint patches of X-ray emission surrounding the radio 
lobes in Hercules A and the bright spot close to the head of the 
western lobe qualitatively fit into this scenario. 

The origin of the jet-like feature in the residual X-ray image 
is unclear One might speculate that it is due to synchrotron-self 
Compton emission from the radio jet. The strength of this 'X- 
ray jet' and the absence of any counter jet would imply that 
Hercules A is not oriented in the plane of the sky. This obvi- 
ously is in contradiction to the existence of a double radio jet 
and the overall symmetry of the radio source (Dreher & Feigel- 
son 1984). 

Another interesting feature to note is the close association 
of the brightest knot in the eastern radio jet and the strongly 
enhanced X-ray emission to the south (and, at a much lower 
brightness level, to the north). It almost seems as if the radio jet 
is squeezed in and probably deflected by local enhancements of 
the surrounding gas density. But again an unambiguous inter- 
pretation is difficult in our case, because the maximum of the 
X-ray emission is close to a rather bright stellar object on the 
optical plate, presumably a foreground star, which might con- 
tribute significantly to the X-ray flux. 

Since X-rays from radio lobes due to inverse Compton 
scattered cosmic microwave background (CMB) photons are 
mandatory, we also estimated the amount of X-ray emission 
expected from this effect. Applying the formalism described in 
Feigelson et al. (1995), we first determined the magnetic field 
within the radio lobes from equipartition arguments. As a con- 
servative lower limit we get Ri 10/zG. Using this value for the 
magnetic field we estimate that at most a flux of Fo.i-2.4koV ~ 
4 X 10^^^ erg cm^^ s^^ is expected from this process. Since 
this flux level is about a factor of 20 below the background flux 
level in our observation, there is no chance to detect these X- 
rays. 

4. Summary and conclusions 

We investigated the spectral and spatial properties of the X- 
ray emission of the prominent radio galaxy Hercules A using 
ROSAT HRI, PSPC, and ASCA observations. The ASCA data 
clearly favor a thermal plasma model with a temperature of 
about 4.3 keV and abundances « 0.4 solar. A small (w 8 per 
cent) non-thermal contribution of the central AGN is consistent 
with the ASCA data. 

The X-ray source associated with Hercules A is clearly ex- 
tended in the ROSAT HRI observation. In addition, it is slightly 
elongated in the direction of the radio jets and lobes. The outer 
part of the surface brightness profile perpendicular to the radio 
structure is well described by a King-model with /3 = 0.63 and 
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Fig. 5. Radio contours overlayed onto the HRI image after subtraction of the best-fit /?-model. The radio map was obtained with 
the VLA at 5GHz (Dreher & Feigelson 1984). 



a core radius of about 120 kpc. The total mass within 600 kpc is 
8.4 X 10^^ Mq and the gas mass fraction is about 18 per cent. 
A central point source contributing about 8 per cent to the total 
X-ray flux is clearly needed to fit the surface brightness distri- 
bution. We determine a luminosity of io.i-2.4keV = 3.4 x 10^^ 
erg s^i for the AGN and Lo.i-2.4keV = 4.3 x 10^^ erg s~^ for 
the extended emission. 

After subracting the extended X-ray emission, significant 
residuals remain in the X-ray image, which partly coincide with 
the radio emission. They strongly suggest an interaction of the 
radio jet with the intracluster gas, an effect which has already 
been noted previously for Cygnus A (Carilli et al. 1994), NGC 
1275 (Bohringer et al. 1993) and MRC 0625-536 (Otani et al. 
1998). 

Spectral and spatial analysis thus provides clear evidence 
that Hercules A is indeed located in a cluster of galaxies, but 
spectroscopic observations are needed to unambigously con- 
firm this interpretaion. The non-spherical shape of the extended 
X-ray emission indicates that the cluster is not relaxed yet, a 
view which is supported by the identification of a double nu- 
cleus in the central cD galaxy, thus indicating an ongoing merg- 
ing process. 
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